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Abstract 

The vacuum-driven nonperturbative factors Li for quark and gluon Green's func- 
tions are shown to define the nonperturbative dynamics of QGP in the leading 
approximation. EoS obtained recently in the framework of this approach is com- 
pared in detail with known lattice data for /i = including P/T'^, £/T^, 
The basic role in the dynamics at T < 3Tc is played by the factors Lj which are 
approximately equal to the modulus of Polyakov line for quark Lfund and gluon 
Ladj- The properties of Lj are derived from field correlators and compared to lattice 



data, in particular the Casimir scaling property Ladj = {Lfund)^'^'-^'^"''^ follows in 
the Gaussian approximation valid for small vacuum correlation lengths. Resulting 
curves for P/T"^, e/T'^, ^"t^ are in a reasonable agreement with lattice data, the 
remaining difference points out to an effective attraction among QGP constituents. 
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1 Introduction 



Dynamics of Quark Gluon Plasma (QGP) is now of great interest, since nu- 
merous results of heavy ion experiments call for strong and possibly nonper- 
turbative forces between quarks and gluons, which cannot be explained in the 
framework of perturbation theory, see [1] for reviews of recent results and their 
interpretation. 
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Recently one of the authors has proposed a new approach to the study of the 
QGP dynamics [2], where the main emphasis was done on the vacuum fields, 
and the resulting modification of quark and gluon propagators was considered 
as the first and the basic step in the nonperturbative (NP) treatment of QGP, 
called Single Line Approximation (SLA). 

As a result one obtains NP Equation of State (EoS) of QGP in the form of 
free quark and gluon terms multiplied by vacuum induced factors. The latter 
are expressed via the only (nonconfining) colorelectric correlator £)f (a;)[3] and 
happened to be approximately equal to the absolute values of Polyakov loops 
L fundi Ladj for quarks and gluons respectively. 

Thus all vacuum NP dynamics in this approximation is encoded in Lf^nd, ^-nd 
Ladj = {LfundP'^ by Casimir scaling [4]. 

Moreover, the phase diagram was calculated in SLA [5] assuming that the 
phase transition is again vacuum dominated, i.e. a transition from confining 
vacuum with vacuum energy density Sconf — ~%G2{conf) to the nonconfining 
vacuum with s^ec — —^G2{dec). 

The resulting phase curve Tc(/x) in [5] depends on AG2 = G2{conf) — G2{dec) 
and is in good agreement with lattice data for standard values of G2{conf) [6] 
and AG2 ~ 0.35 G2 (conf). 

Thus the SLA is a reasonable starting point with no fitting or model param- 
eters, since Lf^nd can be computed analytically [3,7] or on the lattice [8,9], 
and AG2 is the fundamental parameter of QCD [5] . This picture of the QCD 
phase transition was called in [5] the Vacuum Dominance Model (VDM) orig- 
inally proposed in [10] in a simplified form (sometimes called the Evaporation 
Model). 

In the model the basic element of the NP dynamics of QGP is the quark and 
gluon Polyakov lines, which are connected to each other by Casimir scaling. 
It is the purpose of this paper to study in detail properties of Polyakov lines 
with the help of the Field Correlator Method (FCM) [11] where those can be 
derived from the nonvanishing colorelectric field correlator Df. In particular, 
Df{x) can be derived from the gluelump Green's function, and the latter 
was calculated analytically in [7,12] and on the lattice the gluelump spectrum 
was found in [13] . These properties can be compared to the lattice data both 
at T ^ Tc and T > Tc, and we predict behaviour of L fund, Ladj at T ^ Tc 
which violates Casimir scaling for n/ = 0, since there Lfund = and Ladj — 
exp (— ^) , with m - known gluelump mass. At T ^ Tc Casimir scaling follows 
from the dominance of quadratic (Gaussian) correlator, and we estimate the 
admixture of higher correlators, violating the scaling. 
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At this point we notice, that contribution of bound states of static quark or 
static adjoint charge with gluons in QGP to the lattice defined Fs{oo,T) and 
consequently to L fundi Ladj would violate Casimir scaling, and the accurate 
observation of Casimir scaling in [4] thus poses some limits on those bound 
states. 

As a result we fix the form of Lfund,Ladj based on our analytic and lattice 
calculations and enter with those to compute EoS, i.e. P{T),s{T) and their 
derivatives. Results of these computations are compared with numerous lat- 
tice data and shown to agree reasonably well within the accuracy of lattice 
simulations. 

The paper is organized as follows. In section 2 basic thermodynamic equations 
for QGP are derived, and natural appearance of Polyakov loops in EoS is 
derived. In doing so an economic expression for the gluon pressure is first 
obtained, while that of quarks is taken from [2]. 

In section 3 the expressions for Lfund-,Ladj are derived in terms of field corre- 
lator and finally in terms of the gluelump Green's function and properties 
of L fundi Ladj both bclow and above Tc are discussed in detail in comparison 
with lattice data. 

In section 4 EoS, P(T) for Uf = 2,2 + 1,3, e(T) and nonideality are 
calculated using the formulas of section 2 and compared to the lattice data. 

Section 5 is devoted to the discussion of results and conclusions. 



2 Derivation of EoS for quark-gluon plasma 

Our derivation below is based on the formalism suggested in [2], where the 
Background Perturbation Theory (BPTh) is exploited, originally worked out 
in [14] and developed in connection to the FCM in [15] for T = and in [16] 
for T > 0. Correspondingly one splits the gluonic field ^4^ into background 
part and valence gluon part a^, as 




Ba + a 



(1) 



and writes the partition function Z{B, T) as 



Z{B,T) ^ N [ D(f)ex.p{- f dt f d^xLtot{x,t)) 



(2) 
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where denotes all set of fields a^, V', and ghost fields. In the lowest order 
in gUfj, one obtains the result in the so-called Single Line Approximation 

Z{B,T) = N,[det{G-')]-'/'det{-Dl{B))[det{ml - b\B - (3) 

where A^i is normalization constant, Dx{B) = dx — igBx, — D\5^y + 
2igFni,. In what follows we put //^ = 0, and consider the case //g 7^ in a 
subsequent paper [17]. 

The thermodynamic potential F{T) is connected to Z{B,T) in a standard 
way 

F{T)^-T]n{Z{B,T))B (4) 

where the subscript B in (Z) b implies avaraging over all background fields. As 
a result F{T) in SLA is a sum of ghion and quark degrees of freedom separately, 
F(T)sLA = Fq(T) + FgiiT). In what follows we omit the subscript SLA, since 
all results (except for corrections to Polyakov lines in the next section) will be 
valid in this approximation. Using the Fock-Feynman-Schwinger (FFS) path 
integral formalism (see [18] for reviews) one has a convenient representation 

- /f ^i{s)d'x{Dzr,,e-^ (^ir($^(x, x))b - (ird(x, x))b ) (5) 

Here 

$f (x, y) = PpP exp {ig j'^ B^dz^^ exp {2ig £ F(^(r))dr) (6) 

and $(a;, y) is the same as in (6) without the last factor, ^(s) in the regularizing 
factor, for details see [2,16]. A similar representation for quarks and antiquarks 
looks hke [16] 

1 1 d s o 

7j^F,{T) = --tr -d'xas)iDz):^e-^-^-'W^{C„) (7) 

Note that in (5), (7) is present the "winding path measure", introduced in [16], 
e.g. for quarks 

00 00 d% ( / M 

n E(-l)"7^ exp UpJY: Urn) -ix- y), - n(36,, (8) 

n=l n=0 V / I. \m=l / ) 

And the same for gluons, {Dz)'^y but without the (—1)" factor. At this point 
we are posing to contemplate the structure of our result (5), (7) and recognize 
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that it is a sum of individual quark of individual quark or gluon lines (Green's 
functions in background) over paths from {lc,0) to {lc,nj3). 

It is clear that for T < T^. this contribution should vanish because of confine- 
ment, and one should look into the representation containing gauge invariant 
Green's functions. These come from white systems, e.g. for singlet {gg) or (qq) 
and the corresponding partition function has the form 

Z^^^K^i, ^) = / dr,dr2{trW{Ci^\ )) (9) 

where dTi arc phase space factors. Note the coinciding indices in W{C^\ C^^), 
which denotes the closed Wilson loop (with possible insertions of F^,^ and 
o'pxFpx for quarks) starting at points (xi,0), (^2,0) (connected by a parallel 
transporter) and ending at points (xl, n(3), (a^, n(3) (again connected). Now, as 
was shown in [2] , in the deconfined phase the pair partition function factorizes 
in the leading approximation of (ga^)'^, while the color-electric correlator 
yields nonzero contribution to each quark or gluon in the form of Polyakov 
lines. The derivation is shortly as follows (see [2] for details). 

The Wilson loop in (9) can be calculated in terms of field correlators using 
cluster expansion theorem [19] 



^Jr{W{Ci'\Ci'^)) = i-tr(Pexp(i^ L^^^''^^^ = 



y- 12 -jr J dSf,,^^{ui) . . . ds^^^^{uk)D^^^^_„^^^^{ui, Uk)j (10) 

where C„ is the total closed loop, containing C^\C^^ and parallel transporters 
from xi to X2 and back, Sn is a surface inside Cn, while the field correlators 
are defined as follows, e.g. for A; = 2 (Gaussian approximation) one has 

{ui, U2) = — tr(F^i^i$(tii, U2)F^^^^^U2, ui)) (11) 

In what follows we concentrate on color-electric correlators, which can be 
written in terms of two scalar functions D^{w),Df{w) (for contribution of 
other color-magnetic correlators see [2,20]. Note that latter do not produce 
factorized contribution, but can support weakly bound states with angular 
momentum L > 0) 

DuMi^^y) = ^{trEi{x)^x,y)Ek{y)^y,x)) = 

S.^iD'^ + Df + ul^) + u,Uk^ (12) 
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where = D^{u), Df = Df{u), u^x-y. 

We now take into account according to [7,12] that correlation lengths A"^ and 
Af, defined from asymptotics D^{u) ~ exp (— |m|/A^), Df ~ exp (— |M|/Af ), 
are small, A^, Af < 0.2 fm. Indeed from the gluelump correlators [12,13] it 
follows that A^ ^ 0.08 fm, Af 0.16 fm. Then for temperatures T < 
1/Af the n dependence appears explicitly [2] and one can write 

^Jr{W{Cl^\ C^^)) = exp (-^(2) - w^^ - . . .) (13) 

where the Gaussian contribution is expressed via , Df 

w^^^ = n/? {VD{r, T) + V^{r, T)) (14) 
and we have defined [2,3] 

VD{r, T)^2j^ dv{l - vT) j\r - Od^D'^{y/e + ^') (15) 

V, (r, T) = 2 du(l - vT) ^d^Dfi^^ + (16) 
Here r = \xi — X2\. 

Now it is clear, that in the confined regime, when is nonzero, VD{r,T) 
grows linearly with r, and factorization of Z^^^C^p ,C^^) is impossible - gg 
and qq propagate as hadrons. For T > T^, however, = and 

Vi{r,T)^V,ioo,T) + v{r,T) (17) 

where v{r,T)\j.^^ = and contains both perturbative and NP contributions. 
In [9] it was shown that v{r,T) is able to support weakly bound states of 
heavy quark and antiquark (QQ), as well as (gg) and Qg systems. This is in 
agreement with lattice data [21]. In the SLA approximation we neglect in the 
first step the effect of f (r, T) and keep only Vi(oo, T). As will be shown below 
this latter contribution explains EoS of QGP with good accuracy. Then the 
gauge invariant quark- antiquark Green's function factorizes into a product of 
one-body terms, each obtaining a factor 

i?L^exp(-«%I^) = (LPy" (18) 

For the gluon gg system one obtains in addition in the exponent the Casimir 
factor = f , which follows from (10)-(12), when all fields are in the 

adjoint representation, 

i2 = exp(-.?li^) (19) 
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In the next section we study these factors in more detail and estabhsh their 
relation to the Polyakov loop factors measured on the lattice. We end this 
section with the discussion of higher correlators in (10). 

Keeping for smooth surfaces only even power correlators (see discussion in 
[22]) one can estimate the contribution of the k = 4 correlator compared to 
the Gaussian one in the exponent of (10) as giving additional factor 

7,; (4) n o npV^ 

- - " - (gm^^'^r ^ M^E? - 'r^Z. < 0-1 (20) 



Wh 



^2) v^- / V- / (2 GeV)2 



where in confinement phase we have estimated (gF)'^ from the string tension 
aE = iJd^xD^ix) ^ XgFY {X^ ^ 0.2 GeVl Note that estimate of {gFf 
from the gluonic condensate yields t] order of magnitude smaller. The estimate 
(20) gives a reasonable explanation of the good accuracy of Casimir scaling in 
the confined phase (see [23] for discussion). In the case of deconfinement, when 
(Te = and gluonic condensate is roughly twice as small (up to T ^ 1.5Tc [8]) 
the "Casimir expansion parameter" t] should be even smaller, since Af does 
not change significantly [8], while , \^ stay constant in this region] ^ I 

We obtain doing path integrals as explained in [2], 

P,-^^ r-Ve-4f^L(") (21) 

^« = "S L I'""'"' S (22) 
Finally performing integration over ds one has 

P,' = ""-^^T' t ^LZ, (23) 

" n=l "' 



_ AN.Uf ^ ^ {-I, 

" n=l 

where we have defined Ladj = L^alji Lfund = L^fundy see (18), (19), and 

= ^ r ^e-?^-4f^ = !Mkj!^) (25) 
WT^Jo 2T2 T ' ^ ' 

These equations and another, integral form instead of the infinite sum, will be 
used in section 4 to compare with lattice data. 



\n+l 



^ Note that estimate (20) refers to the plane or at least smooth surface, while for a 
crumpled surface higher correlators play important role to make area law with the 
minimal surface. 
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3 Polyakov lines and field correlators 



Below only quadratic (Gaussian) field correlators are considered, basing on the 
Casimir scaling property which these correlators ensure, and being in agree- 
ment with lattice data both for T = [23] and for T > [4]. At T > 

four Gaussian correlators are (x) , Df (x) , (x) , (x) , with a^'^ = 

1 J D^'^ {x)(fx. At T > Tc the correlator and vanish, as was suggested 
in [10] and proved on the lattice [8], and three other correlators are nonzero, 
moreover the spatial string tension ag = (t^ grows with temperature in the 
dimensionally reduced limit [24]. This fact explains also the growth with tem- 
perature of the Debye mass, m/j = I^JWl [25], which is known from lattice data 
[24]. Apart from this quantity, we shall not use below the colormagnetic cor- 
relators, since they do not produce static potentials for interparticle angular 
momentum L = 0. 

Therefore we shall be interested only in color-electric (CE) correlators D^{x) 
(inside the confining phase bounded by the curve Tc{f^)), and Df{x) in the 
whole fjL, T plane. 

It is important at this point to stress that in our approach only gauge invariant 
states |n) are to be considered in the partition function at T > 0, 

Z = ^(n|e-^/^|n) (26) 

n 

as well as in all QCD states at T = 0. This is evident in the confining phase, 
since a colored part of the gauge invariant system is connected by the string 
to other parts. 

With the lack of string in the deconfined phase the necessity of using the 
gauge invariant amplitudes is less evident, except for worldlines in the spatial 
directions, where colormagnetic confinement with nonzero ag is operating. 

Nevertheless our use of gauge invariant amplitude, which factorizes at large 
interparticle distances in the deconfined phase, leads to the explicit prediction 
of EoS with modulus of phase factors, which approximately equal to modulus 
of Polyakov lines. 

Below we shall use, as in [2,5], the gauge invariant states, |n) at all fi, T and we 
shall express the interparticle dynamics in terms of gauge invariant quantities, 
like pair or triple static potentials. The large distance limit of these potentials 
yields one-particle characteristics - the self-energy parts of quarks, antiquarks, 
gluons etc. One can use those to study thermodynamics of QGP in the one- 
particle, or Single Line Approximation (SLA) [2]. It is rewarding, that the 
field correlator method is a natural instrument in describing this deconfined 
dynamics, since in absence of the correlator Df has the form of the full 
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derivative and produces gauge invariant one-particle pieces - self-energy parts 
- automatically (in addition to interparticle interaction decreasing at large 
distances) . 

The gauge invariant states \n) . {n\ formed with the help of parallel transporters 
(Schwinger lines) ^{x,y) = P exp{ig Jy Afj_dz^), create, as shown in [2], Wilson 
loops W{C) for qq, qqq, or else (qqqq) systems. From the latter, as shown in 
[2,3], one obtains static potentials. When treating colored systems like (qq), 
the latter is taken as a part of gauge invariant system (qqqq), and the pairs {qq) 
and qq) are separated at large distance where potential V{qq, qq) is neglected. 

We start with the color singlet qq system and write contributions of D^, Df 
at nonzero T — 1/ 13 to the static potentials [3] 

V,{r, T) = dv{l - vT) ^d^Df{^e + ^') (27) 

Voir^T) = 2iyiy{l - vT) {r - E,)diD^ {^e + ^^) (28) 

It is important that Vi, Vd give the contribution to the modulus of Polyakov 
loops, namely [3] 

ij::L = exp(-iiP]±^), lS = (0"^ (29) 

where Vi{T) = Vi(oo,T), Vd = VD{r*,T) and r* is an average distance 
between the heavy quark line and light antiquark (for n/ > 0) , or "heavy 
gluon hue" and a gluon for Ladj- The Casimir scahng relation (29) predicted 
in [3] is in good agreement with lattice data [4], as well as vanishing of Lfund 
for T ^ Tc, n/ = and the strong drop of Ladj for T ^ Tc- Indeed, for T ^T^. 
and Uf = one has r* — > oo and Vd — > cxo, explaining the vanishing of L^f2nd- 
For L^^j in this region one can take into account the kinetic energy of the 
gluon in the system adjoint source plus gluon in a gluelump. This yields an 
estimate Ladj{T ^ Tc) = exp (— ^^), where nigip was computed in [12,13] to 
be K 1 GeV. 

In [3] it was mentioned, that Polyakov lines measured repeatedly on the lattice, 
are expressed through the (singlet) free energy of QQ system at large distances 
Fqq{oo, T) in the same way as in (29), i.e 

Lfld = exp ( 1 (30) 

and actually the difference between Lj-^"* and Lj-^'* was not taken into account 
in [3]. This difference can be easily seen in the standard representation of 
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where n{QQ) denote all excited and bound states where QQ participate, and 
V^^{r,T) is the energy term of such state n when distance between static 
charges Q and Q is equal to r. It is clear that L^^'* coincides with L^p when 
all states n except for the ground state n = are neglected. In this case 
V'o*^*^(r, T) coincides with Vi(r, T), and hence with ^^^(r, T). Note at this 
point, that Vi(r, T) in (27) does not depend on T in the limit when the vacuum 
correlation length X{Df{x) ~ e"''/^) tends to zero, T < ^ « 1 GeV. 

In the general case all states n{QQ) contribute and therefore (c^ > 0) one has 
inequality 

V,{t,T) ^ FIq{t,T) (32) 

To define Vi and Lfund properly, one should separate perturbative and NP 
parts and renormalize Vi to get rid of perimeter divergences. 

The separation in Df{x) can be seen at small x [7] 



Df{x) = 

I l + 0(^ln^ + ^ + . . .1 = --(x) + Dt'\x) (33) 
and at large x, Di^^\x) is [7]. 



Di'^\x) ^ A^^^, A, = C{f )a,2Moaadj (34) 
V x'^ 

where Mq is the lowest gluelump mass [12,13], Mq ~ 1 GeV. 

The corresponding separation of Vi{r,T) is done in [3,9] as follows 

V,{r, T) = Vr\r, T) + v/"^) (r, T) + v/'*"^ (a) (35) 

where 

V^'''\r, T) = _^iQ^e-mo'^(i + o{rT)) (36) 
is as in (27) with Df Dl^{x), 

yidiv),. ^ 2C{f)as n ^ o(rina)') (37) 
n \a J 
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Fig. 1. Shown on the figure are curves of Ladj (blue dashed) and Lfy_nd (red 
dashed) compared to the ones taken from [4]. In the T < Tc region the 
M(q!<j = 0.195) = 0.982 GeV gluelump mass was used. In the deconfinement re- 
gion the fit (39) was used with = 270 MeV for Lfund and the Casimir scaled 
value for Ladj- 

Here rriu = mD{T) ^ 'info's is the np Debye mass [25], and a is the lattice 
cut-off. 

The renormalization procedure suggested in [3] amounts to discarding vf™^{(i)^ 
and this is in agreement with the lattice renormalization used in [26], where 
F^Q{r,T) was adjusted to the form Vf^''*(r, T) at small r and T. Note, that 
V{^^{r,T) ~ O(r^) in this region and the procedure indeed allows to eliminate 
the constant term Vf^'"{a). 

We start with the one-particle limit of Vi(r, T), and the corresponding contri- 



bution to L 



(v) 

fund' 



According to the discussion above, one defines the renormalized Polyakov loop 
as in (29), (35) with Vi(T) = V"^(oo,T) and we shall neglect the difference 
between Lj and L^^^ (important at large T, where Lj^^^ > 1, while always 



From (34) one has (at T ^ Tc 



fl 



Vr'(oo.T) = ^ 



1 ^ f 1 

1 1 — e T 

Mo 



(38) 



so that \//"^V,Tc) ^ ^^^^^^ ~ 0.5 GeV for Mq ^ 1 GeV [12,13]. 
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The same type of estimate one obtains from lattice data [27] where at T > Tc 
one can parametrize the data as follows 

^qq(oo,T) ^ , , ,^qq(oo,T,) ^ 0.5 GeV (39) 

1.35 (^.J - 1 



Thus one can say that quarks (and antiquarks) have selfenergy parts Kq{T) ~ 
Kq{T) = lVi{T) ^ IF^q{oo,T) ^ 0.25 GeV at T fa Tc. 

To illustrate our discussion of Vd, V\ and Lf^nd, Ladj we show in Fig.l our 
curves for L fundi Ladj computed from (29) with Vi(cxd,T) = Fqq{oo,T) taken 
from (16) for T > while L^dj = exp (-^) for T ^ T^. Our dashed curves 
are plotted in Fig.l in comparison to lattice data from [4]. 

For gluons one has instead Kg{T) = ^Kg ^ 0.56 GeV. Let us turn now to the 
r-dependence of interaction. The perturbativc part has a standard screened 
Coulomb behaviour (36), while the NP part vanishes at small r; 

yr(r, T) ~ const • r^ r^O (40) 
From (27), (34) one has as in [3] 

T) - ^MMor)Mor + 0{^) ^ V^/"^ V, T) + ^(r, T) (41) 

Hence the NP interaction in the white system QQ changes from V]"^(oo, T) ^ 
0.5 GeV at large r to zero at small r. The same (multiplied by |) is true for 
the white gg system. 

We end up this section by discussion of the role of excited states in definition 
of F^Q and possible violation of Casimir scaling for Lf^nd, Ladj- It is clear that 
in Fqq for Uf = the only possible excited states consist of gluons {Qg){Qg)', 
{Qgg){Qgg) etc. As it was shown in [9], the weakly bound states {Qg) indeed 
are supported by Vi(r, T), and neglecting the small binding energy the total 
energy of these states is roughly the sum of selfenergy parts kq and Kg 

Eq9 ~ \yA^.T) + ^1/1(00, T) ^ 0.8 GeV(r « T,) (42) 

This should be compared to the possible bound state of an adjoint static source 
G plus gluon, which in the weakly binding limit can be written as 

~ 2 • Vi(oo, T) ^ 1.1 GeV(r ^ (43) 
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In addition multiplicities of states (42) and (43) are different, which leads 
to different predictions for corrections to F^q and Fqq, not connected by 
Casimir scaling, in contrast to the main (ground state) term V^q = Vi(oo,T) 
|Vi(cx),T). Therefore one expects violation of Casimir scaling by 



and Vqq 



gluon- induced bound states in Lfund and Ladj, and high accuracy of lattice 
data [4] indicates then a small role of such bound states. 



4 A comparison to the lattice data 
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Fig. 2. Pressure ^ as function of temperature T . Shown on the left figure is a 
comparison of the analytical calculus (48), (49) (dashed lines) with the lattice results 
(bold lines) [28] for the case = 0,2,3. Shown on the right figure is the case of 
Uj = 2 + 1. Green dashed line is the analytical calculation (48), (49) compared to 
the lattice one from [29]. 
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Fig. 3. Pressure ^ as function of temperature T. The case of nj = 2 + 1 (left figure) 
and Uf = 3 (right figure) (48), (49). Lattice results were taken from [30]. 



In this section we shall exploit the reduced pressure p = which for /z > 
can be written as: 
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(45) 



with 99^") given in (25) 



2 2 
n mt 



in) ^ l^Hlj^ I !^ 



(46) 



Both sums can be written in a more convenient way. Using the representation 

of i^2, 



(n) 



where u = rUq/T, one haS: 



6 7o TFT 



(47) 



3 7r2 



( S - - j + (^ag + - 



(4J 



- 1 /-"c z^dz 



37r2 Jo e^+'^s' - 1 



with ttg = Vi(T)/2T, ttgi = |aq and 



(49) 



/■OO ^4 

In the paper we consider the case of n = and characteristic temperature 
region of T ^ [T^ = 170 -7- 270 MeV) where quark masses do not affect the 
thermodynamical functions appreciably. This is due to the fast convergence 
of the sum over n at large n ensured by factors L" (L < 1) while ip^^^ ~ 1 
for n 1. Characteristically, Lp^^\mq = 0) = 1, and for = 0.4T one has 
Lp^^^ = 0.96, (f^g^^ = 0.03. Therefore one can with a good accuracy neglect 
masses in (48), (49): 

2nf f°° z^dz 

P<i=2 z+a 11 ^1 

TT^ JO e^+"9 + 1 

8 z^dz 
^^'=3^/0 ^''^ 

Eqs. (51), (52) are compared with lattice pressure data in Fig.2 for rij = 2 + 1 
(left) and Uf = 3 (right figure). In Fig. 3 are shown our calculated curves for 
the cases nj = 2 + 1 (left part) and Uf = 3 (right part), which are compared 
with lattice data from [30]. 



The form (48) was independently obtained by N.O. Agasian (to be published). 



14 



To simplify further one can use for /i = instead of (51), (52) the first terms 
of expansion in (44), (45), namefy: 



12n 



f 



71 

16 

TT 



■'fund 



(53) 



P9 = :;^^adj (54) 

Another useful quantities to compare with lattice data are the internal energy 
density and the "nonideality" of the QGP: 



(55) 




This work: nf=2+l: d=0.1ms and ^=0.2m^ 



O Nj=6, m^,d=0.1m5 
O Nt=6, m,i,d=0.2m, 
A Nt=4. m^,d=0.1mi 
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pA, Ni=4, nf=3 1- 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 



Fig. 4. Energy density ^ as function of temperature T. The case nj = 2 + 1 with 
Tnu^ = 0.1ms and rriu^d = 0.2ms (green dashed curve) (55) is compared to lattice 
data from [29](left fig.). The case nj = 2 + 1 with m^^^ = OAT, nig = T (red dashed 
curve) and nj = 3 with m^ = 0.4T (blue dashed curve) (55) are compared to lattice 
data from [28] (right fig.). 

Using (48), (49) one has 



dz 



(56) 



-(0) 



and the "nonideality" of the QGP: 



e^+"fl' + 1 



e-3P 

rp4 



dT 



(57) 



(5J 



In the simple approximation (53), (54) one has 

/(T^) = ^'^^f rp dLfund ^ IQrpdLadj 



dT 



7r2 



(59) 
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m,= 2MeV, md=6MeV, ms=100MeV 
T — ] — I — I — I — I — j — I — I — r^i — I — I — ! — T — I — I — I — I — r — I — 1^ 




Fig. 5. Energy density ^ as function of temperature T. The curve for n/ = 3 with 
ruu = 2MeV,mrf = 6MeV,ms = lOOMeV (green dashed) (55) is compared to lattice 
data from [31]. 

We compare our calculations for ^ in Fig.4 and 5 with three different lat- 
tice data: [28], [29], [31]. In Fig.6 we demonstrate our I(T) computed from 
(58), (51), (52) with lattice data of 2+1 flavor from [31] (left curve) and from 
[29] (right curve). 



mu=2MeV, rnd=6MeV, m,=100MeV m„,j-0.1m, and m,,j-0.;m, 




1 1.5 2 2.5 3 Ti""^^! 



T/T„ 

Fig. 6. "Nonideahty" of QGP (e-3p)/rl Shown are the curves for (left fig.) nj = 3 
with rriu = 2 MeV, = 6 MeV, = 100 MeV (green dashed line) compared to 
[31] and (right fig.) for n/ = 2 + 1 with niu^ = 0.1ms and mu^d = 0.2ms compared 
to [29]. Analytical calculations are done using (48), (49), (55). 
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At this point it is instructive to estimate the contribution of qq. gg interac- 
tions to the pressure. Writing the virial coefficient in the form Pj = Pj^\l + 

-^Bj{T) + ...), where P^ = Pg, Pgi in SLA, Eqs. (21), (22), with 

Bj{T) = ^ /(I - e^^(^'^)/^)cil-, J = fund,adj (60) 

and taking for qq and gg interaction term U fund and Uadj respectively at large 
T as Uj{r,T) — Tuj{rT), one obtains a corrected pressure 

P = Pi'\^ - c,) + P^'!\l - c,i) (61) 

where Cgi ^ ^ /o°^ p2^p(el""*('')l-l), c, ^ /o°^ p2^p(el"/'""'('')l-l). Note that 
qq and gg interaction in the singlet color state is attractive, so that \Uj\ — —Uj. 
The dependence on rT in uj occurs at large T, in the dimensionally reduced 
regime, when dynamical dimensional quantity is the spatial string tension 
an = const • T^, and the Debye mass moiT) = I^/oh = const • T. 

Thus one expects that 1) the corrected pressure is smaller than the SLA 
predicts, 2) the large T behavior of P{T) is below the Stcfan-Boltzmann values 
(modulo logarithmic factors). Both features are clearly seen in the Fig. 2, 3,4, 5. 



5 Discussion of results. Conclusions. 

We have shown in section 2, following [2], that EoS in the zeroth approxi- 
mation is represented by free quark and gluon lines augmented by the factors 
Lfund for quarks and Ladj for gluons. These factors have been derived from the 
Gausssian color-electric correlators D^{x), Df{x), and the latter in its turn 
can be computed analytically from the gluelump Green's function, or directly 
on the lattice [8,9]. This representation of Ladj and Lf^nd allows to express 
Lj, j — fund,adj in terms of the NP static potential Vi{r,T) at r = oo, 
and compare the latter with the singlet free energy FgQ{r,T). It was argued 
that Vi and Fqq differ due to presence of excited Qg^ states in Fqq, and can 
be taken equal in the first approximation. This leads to the identification of 
Lj with the modulus of corresponding Polyakov lines. In the Gaussian ap- 
proximation for Vi one then automatically obtains the Casimir scaling for Lj: 

Ladj = (-^/und) '^2(/""<*) which is observed on the lattice with good accuracy [4]. 
Corrections are found to be of two types: 1) contribution of higher correlators 
to Vi and Lj yields less then 10% (20) and can be neglected 2) contribution of 
excited states yields corrections not connected by Casimir scaling and there- 
fore high accuracy of data [4] imposes a stringent limit on the role of excited 
states of the type {Qg"')- For T > our expression (29) automatically predict 
vanishing of Lf^nd for n/ = and behavior of Ladj — exp {—Mq/T) with Mq 
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- lowest gluclump mass ^ 1 GeV. These features are in good agreement with 
the lattice data [4], and are shown in Fig.l. 

For EoS using formulas of section 2 and treatment in [2] we have given two 
types of expressions for the pressure P: 1) as a sum over winding n (Matsubara 
frequencies) in (44), (45) and equivalent forms as integrals over "momentum" 
z in (48), (49). It was argued that for = and not large T, ^ T < 2Tc 
one can use much simpler forms of (53), (54), which are first terms of the sums 
(44),(45). 

In all these forms the only source of non-perturbative dynamics in EoS is 
Polyakov factors Lj, which are defined independently and therefore our EoS is 
the explicit prediction without any model of fitting parameters. Hence check 
of our approach is the check of our basic principle that non-perturbative dy- 
namics enters in the form of vacuum based factors Lj. 

Comparison of our EoS, (44), (45) or (48), (49), is done with several lattice 
groups for each quantity, to have an idea of accuracy of our results and of 
lattice data, and dependence on quark masses. The latter appears very weak 
in EoS, e.g. quark mass of rUq — OAT yields a 4% correction to the zero mass 
result, while on the lattice this dependence is stronger. We compare pressure 
Y4 for nf = 0, 2, 3 and nig = OAT in Fig. 3 (left part). One can see deviation of 
~ 20% of our curves from lattice data [28] for T < 3Tc and the same type of 
agreement ior Uf — 2 + 1 with data from [29]. Typically our curves are higher 
with the fact that the (attractive) interaction between quarks, antiquarks and 
gluons is not taken into account. The first correction (60), (61) treating this 
attraction between qq and gg, has the negative sign, which might improve the 
agreement. The agreement is however better with another set of lattice data 
from [30] done for Uf — 2 + 1, see Fig. 3. Comparing left and right parts of 
the Fig. 3 one can notice, that lattice data [30] are much more sensitive to the 
quark masses, than our prediction. 

Another interesting comparison is for the internal energy e and non-ideality 
/ = ^^i^, given in Fig.4,5 and 6. It is important that both quantities contain 

derivatives T '^^^J^^ and therefore are much more sensitive to the type of non- 
perturbative dynamics, which is present in our approach. The agreement of 
our Eqs. (56), (57) with data from [28] and [29] are shown in Fig. 4 and is of 
the same quality as for the pressure: one has ~ 15% higher theoretical curve 
for T > 1.2Tc, the same one can see in Fig.5 with data from [31]. Note, that 
the quark masses in this case are close to physical ones. Finally, the non- 
ideality is compared to the data from [31] in the left part of Fig. 6 and is in 
good agreement with data [29] (right part of Fig. 6 is less successful, because 
lattice data from [29] and [31] differ strongly). As a whole, it is surprising that 
such simple approach without any parameters (actually primitive formulas 
(53), (54) already have sufficient accuracy within our approximation) yields a 
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reasonable agreement with lattice data for P{T), e{T) and I{T). If one adds to 
that a good agreement of our phase curve Tc(/i) in [5] with majority of lattice 
data, the possible conclusion might be, that our zeroth approximation to the 
non-perturbative vacuum fields - taking non-perturbative contribution in the 
form of Lj - is a viable spproach to the dynamics of QGP. The next step is an 
account of possible perturbative and non-perturbative interactions between 
quarks, antiquarks and gluons, which is partly done in [3,9] for color-electric 
fields {Vi{r,T)) and in [20] for color-magnetic ones. The exact contribution of 
these effects to the EoS is not yet done and should be an important next step. 
The strong interaction in (qq) and (gg) systems discovered in [20] might give 
further support for the idea of strong quark-gluon plasma - sQGP. 

The authors are indebted for useful discussions to members of ITEP and FIAN 
physical seminars. We are also grateful to S.N.Fedorov for providing us with 
a useful program " GetData" . The financial support of the RFFI grant 06-02- 
17012 and NSh-843.2006.2 is acknowledged. This work was supported by the 
Federal Agency for Atomic Energy of Russian Federation. 
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